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SUMMARY 

The  reflection  of  a pulsed  lasejr  beam  from  a sea  surface 
depends  on  the  diameter  of  the  beam  and  on  the  wave -number 
spectrum  of  the  water  surface.  General  expressions  are 
derived  for  the  influence  on  the  reflection  of  the  mean 
height  and  slope  of  waves  in  the  illtiminated  area  as  well  as 
the  variation  of  the  height  and  slope  within  this  area.  The 
results  are  applied  to  applications  of  airborne  lasers  for 
wave  height  profiling,  altimetry,  and  water  depth  sounding. 
The  optimum  beam  diameter  depends  on  both  the  application  and 
the  wave-number  spectrum  of  the  wjiter  surface. 


Approved  for  Public  Release. 


POSTAL  ADDRESS:  Chief  Superintendent,  Electronics  Research  Laboratory, 

Box  2151,  G.P.O. , Adelaide,  South  Australia,  5001. 


UNCLASSIFIED 


ERL-0015-TR 


TABLE  OF  CONTENTS 

Page  No. 

1.  INTRODUCTION  1 

2.  DESCRIPTION  OF  THE  SEA  SURFACE  1-3 

2.1  Wave-height  spectrum  1-2 

2.2  Wave-slope  spectrum  2 

2.3  Equilibriiom  spectra  2 

2.4  Wave  statistics  3 

3.  PROPERTIES  OF  A LIMITED  AREA  OF  THE  SEA  SURFACE  3-8 

3.1  Variance  o£  mean  height  3-5 

3.2  Variance  of  residual  height  5-6 

3.3  Variance  of  relative  central  height  6-7 

3.4  Variance  of  mean  slope  7-8 

3.5  Variance  of  residual  slope  8 

4.  LASER  BEAM  REFLECTION  FROM  THE  SEA  SURFACE  8-15 

4.1  Average  amplitude  of  reflected  pulse  8-10 

4.2  Variability  in  amplitude  of  reflected  pulse  11  - 14 

4.3  Water  height  resolution  14  - 15 

5.  DISCUSSION  OF  APPLICATIONS  15  - 16 

6.  CONCLUSION  16 

7.  ACKNOWLEDGEMENT  16 

REFERENCES  17 

LIST  OF  FIGURES 

1.  Dependence  of  the  relative  variance  of  the  mean  height  on  the 
dimensionless  radius  of  the  spot. 

2.  Dependence  of  the  relative  variance  of  the  residual  height  on 
the  dimensionless  radius  of  the  spot. 

3.  Dependence  of  the  relative  variance  of  the  central  residual 
height  on  the  dimensionless  radius  of  the  spot. 

4.  Dependence  of  the  relative  variance  of  the  me£m  slope  on  the 
dimensionless  radius  of  the  spot. 

5.  Dependence  of  the  relative  variance  of  the  residual  slope  on 

the  dimensionless  radius  of  the  spot. 

6.  Schematic  diagram  of  laser  beam  reflection  from  the  sea  surface. 

7.  Geometry  of  tilted  laser  beam  and  receiver  field  of  view. 


ERL-0015-TR 


8.  Dependence  of  the  relative  signal  return  on  roll  angle  and  beam 
divergence. 

9.  Dependence  of  the  relative  threshold  signal  return  on  the  rms 
residual  slope  for  a failure  probability  of  0.1%. 

10.  Dependence  of  the  relative  thresjhold  signal  return  on  the  rms 
residual  slope  for  a failure  probability  of  1%. 

11.  Dependence  of  the  relative  threshold  signal  return  on  the  rms 
residual  slope  for  a failure  probability  of  10%. 

12.  Dependence  of  the  relative  threshold  signal  return  on  the  spot 
diameter  for  three  failure  probabilities. 

13.  Dependence  of  the  relative  threshold  signal  return  on  the  spot 
diameter  for  a failure  probability  of  1%  and  a variety  of  sea 
states . 

14.  Range  of  sea  states  for  which  given  spot  diameters  and  relative 
thresholds  achieve  a failure  probability  below  1%. 

15.  Range  of  sea  states  for  which  given  spot  diameters  and  relative 
thresholds  achieve  a failure  probability  below  0.1%. 

16.  Variation  of  water  height  and  slant  height  errors  over  illuminated 
spot. 

17.  Dependence  of  rms  residual  height  error  on  laser  spot  diameter  and 
sea  state. 


1 


ERL -001 5 -TR 


1 . INTRODUCTION 

The  reflection  of  a short  laser  pulse  from  a sea  surface  is  influenced  in 
several  ways  by  the  wave-number  spectrum  of  that  surface. 

The  mean  water  level  within  the  illuminated  area  determines  the  delay 
between  the  transmission  of  the  laser  pulse  and  the  detection  of  the  peak  of 
the  reflected  pulse.  The  variation  of  mean,  water  level  in  a series  of 
measurements  is  the  data  used  by  a wave  height  profiler.  This  same  variation 
is  also  an  error  that  limits  the  accuracy  of  an  altimeter.  Clearly,  these  two 
applications  require  beams  of  different  diameters. 

The  variation  of  water  level  within  the  illuminated  area  determines  the 
temporal  spread  of  the  reflected  pulse.  Since  this  spread  slows  the  pulse 
risetime,  amplitude  fluctuations  introduce  timing  errors  when  threshold 
detection  is  used.  A small  beam  diameter  will  minimize  water  level  variation 
within  the  illuminated  spot,  but  it  also  increases  amplitude  fluctuations  in 
the  reflected  pulse  detected  by  a receiver. 

The  slope  distribution  of  the  sea  surface  affects  the  amplitude  and 
variability  of  the  received  signal.  When  a rippled  water  surface  is 
illuminated  by  a large  diameter  laser  beam,  the  reflection  is  diffuse  and  the 
amplitude  of  the  received  signal  is  relatively  stable.  However,  when  a long 
smooth  swell  is  illuminated  by  a narrow  laser  beam,  the  reflection  is  almost 
specular  and  dances  around  in  all  directions  producing  an  occasional  bright 
flash  at  a fixed  observation  point.  For  reliable  detection  of  the  reflected 
pulse,  a large  beam  diameter  is  therefore  preferable. 

These  and  other  considerations  will  sometimes  lead  to  conflicting 
requirements  for  the  beam  diameter.  The  theoretical  study  of  laser  beam 
reflection  from  an  ocean  surface  by  Swennen(ref . 1)  takes  no  account  of  the 
wave-number  spectrum  of  the  sea  surface.  In  this  report,  general  expressions 
are  derived  for  the  influence  of  wave-number  spectra  on  laser  beam  reflection 
and  the  results  are  applied  to  several  important  applications. 


2.  DESCRIPTION  OF  THE  SEA  SURFACE 


2.1  Wave-height  spectrum 

The  undulation  of  a water  surface  can  be  described  by  the  height  f (r) 
of  the  water  above  a point  r in  the  horizontal  plane  of  mean  water  height. 
By  this  definition,  then. 


f (r)dr  = 0 


The  surface  heights  at  adjacent  points  are  related  by  the  two-dimensional 
spatial  autocorrelation  function 


Z(r) 


lim  (7r(f)"‘ 


?(P)  ?(r+P)  dp  . 


Q».-  p<  Q 


(1) 


It  should  be  noted  that  Z(oj  is  simply  the  variance  of  the  surface  height. 
The  wave-number  spectrum  of  the  water-height  variance  is  the  Fourier 
transform  of  the  autocorrelation  function,  namely. 


^(k)  = C2^) 


Z(r)  exp  (-ik.r)  dr 


(2) 


where  k = 2ir/X  is  the  wave-number  corresponding  to  the  wave  length  X. 
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Using  the  polar  co-ordinates  r = (r,a),  a direction -independent  wave- 
number  spectrum  is  defined  by~ 

XCk)  = / ^(k)  kda  , (3) 

which  differs  by  a factor  k from  that  used  by  Phillips  (ref. 2,  ch.4). 

2.4  Wave-slope  spectrum 

The  slope  of  a water  surface  at  a point  r can  be  defined  by  the 
cartesian  components  ~ 


s = 5 f (r)  and  s = 

X X y 

from  which  it  follows  that  s = s = o.  The 

X y 

function  for  slope  is  defined  by 

= lini  (’TCf)-’  /j  S^(P3 

p<  Q 


spatial  autocorrelation 

s^(r+p)dp  , 


(4) 


(5) 


and  the  corresponding  direction-independent  wave-number  spectrum  is 

denoted  by  S (k) . Assuming  that  orthogonal  slopes  are  statistically 
xy 

independent,  it  will  be  convenient  to  use  the  wave-number  spectrum  of 
the  total  slope  variance  defined  by 


sck) 


2 

XX 


Ck)  * Syy  Ck)  . 


(6) 


2.3  Equilibrium  spectra 

The  wave-number  spectrinn  of  wave  height  takes  a simple  form  when  the 
sea  is  "fully  aroused".  This  occurs  when  the  wind  has  blown  over  a fetch 
and  for  a duration  that  are  sufficient  for  the  sea  to  be  in  equilibrium  with 
the  wind.  This  equilibrium  is  established  when  the  energy  supplied  to  the 
waves  is  lost  again  as  crests  break  and  form  patches  of  foaming  white  caps 
or  as  parasitic  capillary  waves  are  generated.  Phillips  (ref .2,  Sec. 4. 5) 
derives  the  fimctional  forms  of  the  wave  spectra  in  the  equilibrium  range 


by  similarity 

cons iderat ions . 

They  are 

X(k)  = 

Bk~^ , 

ko<  k^  ki  , 

(7) 

and 

2(k)  = k^x(k) 

II 

09 

1 

, ko<  k<  ki 

(8) 

where  k©  and  ki  depend  on  v/ind  and  sea  conditions. 

The  validity  of  these  functional  forms  is  supported  by  measurements  of 
the  related  frequency  spectra.  Data  from  numerous  experiments  has  been 
assembled  by  Phillips  (ref  .2,  fig. 4. 8)  and  this  is  supported  by  more 
recent  data  obtained  by  De  Leonibus  et.  al.(ref.3)  and  Peep(ref.4).  The 
value  for  B = 4.6  x lO"^  assigned  by  Phillips (ref. 2)  and  Wu(ref.5)  is 
derived  from  the  measurements  of  slope  variance  by  Cox  and  Munk(ref.6). 
However,  the  value  B = 6.7  x 10""*  was  derived  by  the  author(ref . 7)  from  a 
different  interpretation  of  the  same  data. 
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2.4  Wave  statistics 

The  statistical  distributions  of  water  heights  and  slopes  are  also 
important  for  a full  description  of  the  sea  surface,  because  they  go 
beyond  the  average  properties  such  as  variances  and  spectra.  Since  the 
shape  of  a sea  surface  results  from  the  superposition  of  many  independent 
waves  generated  by  winds  in  different  regions,  the  central  limit  theorem 
predicts  gaussian  distributions  for  wave  height  and  wave  slope. 

Experimental  measurements  of  wave  height,  analysed  by  Kinsman (ref. 8) 
reveal  a distribution  that  is  approximately  gaussian  but  has  a significant 
skewness.  Observed  surface  slopes,  analysed  by  Cox  and  Munk(ref.6)  also 
show  an  approximately  gaussian  distribution.  However,  they  also  report 
an  appreciable  skewness  in  the  windward  direction,  that  results  in  the  most 
probable  surface  slope  being  not  zero  but  a few  degrees.  These  departures 
from  gaussian  distributions  result  from  non-linear  distortions  of  the  wave 
profiles,  as  well  as  interactions  and  couplings  between  different  wave 
components.  In  the  present  study,  the  first  order  approximation  of  a 
gaussian  distribution  is  asstimed,  namely  the  probability 

P(s^,s^)  = (2na^a^)-^  exp  I 

where  o and  a are  the  standard  deviations  of  slope  in  the  crosswind  and 
c u 

upwind  directions  respectively. 

3.  PROPERTIES  OF  A LIMITED  AREA  OF  THE  SEA  SURFACE 

The  applications  of  an  airborne  laser  described  earlier  all  involve  the 
reflection  of  a beam  from  a limited  area  of  the  sea  surface.  Important 
properties  of  the  reflection  (such  as  time  delay,  time  spread,  direction  and 
amplitude)  depend  on  the  nature  of  the  surface  of  the  illuminated  area. 
Expressions  are  now  derived  for:  the  pulse  to  pulse  variation  of  the  mean 
height  and  mean  slope  over  the  illuminated  area  of  water;  the  variation  of 
the  height  and  slope  within  this  area;  and  the  variation  of  the  difference 
between  the  central  height  and  mean  height. 

3.1  Variance  of  mean  height 

The  mean  water  height  over  a circle  of  radius  R,  centred  on  the  point 
r is  given  by 

fR(r)  = (7rie)-»  ' II  ?(P)dp  . (10) 

Ip-r  l<  R 

This  expression  can  also  be  written  as  the  convolution 
fR(r)  = ?C£)*\(£)  - 

where  the  averaging  function 


has  the  Fourier  transform (ref .9) 


Hj^Cy  = (27r)-2  2J,  (kR)/kR  . 


(13) 


The  autocorrelation  function  for_the  mean  water  height  is  derived 
by  substituting  the  expression  for  ii'  equation  (11)  in  place  of 

? (r)  in  equation  (1) . This  yields 


Zj^(r)  = Z(i:)*hj^(r)*hj^(-r)  , 


(14) 


which  is  a generalization  of  the  theorem  proved  by  Papoulis(ref .10, 
equations  (12-43)).  It  follows  that  the  corresponding  wave-ntimber 
spectrum  is 


^j^(k)  = ‘ 


(15) 


from  which  the  direction-independent  form 


Xj^(k) 


X(k) 


2Ji  (kR)/kR  f 


(16) 


is  derived  by  substituting  Hn(k)  from  equation  (13). 

When  the  sea  is  fully  aroused,  the  wave-number  spectrum  in  equation  (7) 
should  be  substituted  for  X(k)  in  equation  (16).  The  autocorrelation 
function  for  mean  water  height 


ZR(r) 


fki 

k-3 

2Ji (kR)/kR  ^ 

i ko 

— — 

Jo (kr)dk  ^ 


(17) 


is  then  obtained  from  the  inverse  Fourier  transform  in  polar  co-ordinates 


Z(r) 


X(k)Jo  (kr)dk  . 


(18) 


It  follows  that  the  variance  of  the  mean  water  height  relative  to  the 
total  height  variance  is 


Zj^(0)/Z(0) 


2 (ko  -ki 


fki 

k"^ 

iko 

2Ji (kR)/kR 

dk  . 


(19) 


This  ratio  is  shown  in  figure  1 as  a function  of  the  dimensionless  radius 
koR  of  the  illuminated  spot. 

It  can  be  seen  from  the  figure  that,  when  the  water  height  is  averaged 
over  a small  area  (koR»-0),  the  variance  of  the  mean  height  approaches  the 
total  height  variance  (Z„(0)  Z(0)).  This  provides  some  validation  of 

the  calculations,  because  no  averaging  can  be  achieved  with  a vanishingly 
small  spot.  As  the  spot  becomes  larger,  however,  the  averaging  becomes 
more  effective  and  the  variance  of  the  mean  height  is  reduced.  Indeed, 
when  koR  exceeds  about  3.5,  the  mean  height  variance  is  only  about  1%  of 
the  total  height  variance.  The  standard  error  in  the  mean  height,  which 
is  the  square  root  of  the  variance,  is  then  about  10%  of  the  rms  wave 
height.  Further  increases  in  the  size  of  the  spot  produce  only  a small 
additional  reduction  in  the  mean  height  error. 
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It  can  also  be  seen  from  the  figure,  that  the  relative  variance  of  mean 
water  height  depends  primarily  on  the  parameter  koR  and  not  to  such  a 
large  extent  on  the  ratio  ki /ko . The  significance  of  the  former  parameter 
is  more  readily  apparent  when  it  is  expressed  as 

ko  R = ttD/Xo  , 


where  D is  the  diameter  of  the  spot  and  X©  (=  27r/ko)  is  the  wavelength  of 
the  longest  significant  wave  present.  The  reason  for  the  rapid  decrease 
in  the  relative  variance  of  mean  water  height,  when  koR  — is  that  the 
diameter  of  the  spot  is  greater  than  the  longest  wavelength  present  and 
hence  the  height  is  averaged  over  at  least  one  wave  in  every  direction. 

The  insensitivity  of  the  averaging  to  the  ratio  ki /ko  is  due  to  the  strong 
bias  towards  short  wave-numbers  of  the  wave  height  spectrum  (equation  (7)). 
In  physical  terms,  this  exemplifies  the  small  contribution  to  wave  height 
from  ripples. 

In  the  application  of  a laser  (or  radar)  altimeter,  the  influence  of 
swell  on  height  measurements  can  be  largely  eliminated  by  making  the 
diameter  of  the  illuminated  spot  greater  than  the  length  of  the  longest  wave 
likely  to  be  encountered. 

In  the  application  of  a wave  height  profiler,  which  determines  water 
height  with  a beam  of  finite  diameter,  the  measured  variance  will  correspond 
to  the  quantity  Zj^(O) . Therefore,  the  profiler  will  record  over  90%  of  the 

total  variance  Z(0),  for  all  values  of  ki /ko , if  koR  is  less  than  about  0.3 
(from  figure  1).  In  other  words,  the  profiler  will  record  the  height  of 
most  of  the  significant  waves,  if  the  spot  diameter  is  less  than  one  tenth 
of  the  characteristic  wavelength  of  the  sea  being  studied. 

3.2  Variance  of  residual  height 

The  residual  height  of  the  water  surface  at  the  point  r + g,  within  the 
averaging  circle  is  defined  as 

- ?rC£)  . (20) 

The  autocorrelation  function  for  the  residual  water  height  is  derived  by 
substituting  the  above  expression  for  ?pCr,P)  in  place  of  ? (r^  in 

equation  (1).  With  the  aid  of  equations  (11)  and  (14),  this  yields 

= Z(r)+Zj^(r)-Z(r+p^)*hj^(r+p)-Z(r-p)*hj^(r-p)  . (21) 

The  mean  residual  autocorrelation  over  the  averaging  circle,  which  is 
defined  by 

= (rf)-  11  z;cr,p)dp , 

p<  R 


reduces  to 


f Cr)  = Z(r)  -Z  (r)  . (22) 

lx  rx,/  Ix 

This  is  a generalization  of  the  well-known  statistical  relationship,  that 
the  sample  variance  is  the  difference  between  the  true  variance  and  the 
variance  of  the  mean. 
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This  important  result  confirms  the  behaviour  expected  from  qualitative 
arguments.  That  is,  as  the  illuminated  spot  is  made  bigger,  the  variation 
in  mean  height  between  measurements  is  reduced  and  the  variation  in  height 
within  the  spot  is  increased.  Equation  (22)  shows  that  the  reduction  in 
the  former  variance  equals  the  increase  in  the  latter. 

The  ratio  of  the  residual  variance  of  water  height  averaged  over  the 
illuminated  area  to  the  total  water  height  variance,  namely 

fj^(0)/Z(0)  , 

is  shown  in  figure  2 as  a function  of  the  dimensionless  radius  (koR)  of  the 
illuminated  spot.  It  can  be  seen  from  the  figure  that,  when  the  spot  is 
large  enough,  the  average  residual  height  variance  approaches  the 

total  height  variance  Z(0),  as  would  be  expected.  As  the  spot  becomes 
smaller,  the  residual  variation  within  the  illuminated  area  is  reduced. 

Thus,  when  koR  is  less  than  about  0.1,  the  residual  variance  is  only  about 
1%  of  the  total  variance  for  all  values  of  ki  /ko . 

The  broadening  of  a fast  laser  pulse  reflected  from  the  sea  surface,  due 
to  variation  in  the  water  height  over  the  illuminated  area,  can  be  predicted 
with  the  aid  of  figure  2.  This  broadening  can  introduce  timing  errors  in  a 
wave  height  profiler  or  water  depth  sounder  when  threshold  detection  of  the 
reflected  pulse  is  employed.  In  order  to  maintain  the  resulting  height 
error  below  10%  of  the  rms  wave  height,  figure  2 shows  that  the  spot 
diameter  must  be  less  than  about  3%  of  the  longest  significant  wave 
present  (i.e.  koR<0.1). 

3.3  Variance  of  relative  central  height 

Another  important  property  is  the  variance  of  the  difference  between  the 
water  height  at  the  centre  of  the  illuminated  spot  and  the  mean  height. 

The  autocorrelation  function  corresponding  to  this  case  is  given  by  equation 
(21)  when  p = 0,  that  is 

Z'j^Cr,0)  = Z(r)  + Zj^(r)  - 2Z(r)*hj^(r)  , 


which  yields  the  wave-number  spectriom 


Xr(1c,0) 


X(k) 


1 - 27rHj^(k)  ^ 


When  the  sea  is  fully  aroused,  a derivation  similar  to  that  in  section  3.1 
shows  that 


z'j^C0,0)/Z(0) 


2(kfl-2-ki-")-^ 


fki  I 

— — 

iko  ^ 

l-2Ji  (kR)/kR 

‘dk 


(23) 


This  expression  for  the  variance  of  the  relative  water  height  at  the 
centre  of  the  illuminated  spot  is  depicted  in  figure  3.  As  expected,  the 
variance  of  the  central  point  approaches  the  total  variance  as  the  size  of 
the  spot  increases.  Comparison  with  figure  2 shows  that  the  central 
variance  is  less  than  the  average  variance  over  the  illuminated  spot.  The 
reason  can  most  easily  be  perceived  by  considering  a plane  as  a first 
approximation  to  the  shape  of  the  water  surface  within  the  spot.  The  mean 
residual  variance  over  the  whole  spot  is  then  a measure  of  the  slope  of  the 
plane,  whereas  the  central  residual  variance  vanishes  because  the  centre  lies 
in  the  plane.  It  follows  that  the  residual  variance  of  the  central  point  is 
second  order  and  depends  on  the  curvature  of  the  water  surface. 
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Figure  3 is  relevant  to  the  application  o£  a laser  depth  sounder,  which 
may  receive  a surface  reflection  from  a large  spot  but  have  a narrow 
concentric  beam  penetrating  the  water  surface.  In  this  case,  the  residual 
height  error  of  the  central  beam  is  less  than  (say)  10%  of  the  rms  wave 
height  when 


z'j^C0.0)/Z(0)  < 0.01  . 

Figure  3 shows  this  condition  to  be  satisfied  when  koR  ^ 0.3,  that  is,  when 

the  diameter  of  the  spot  used  for  the  surface  reflection  is  less  than  10% 
of  the  longest  significant  wavelength  present.  Comparison  with  figure  1 
shows  that,  a requirement  for  small  variance  of  the  mean  water  height  over 
a spot  conflicts  with  a requirement  for  small  residual  variance  at  the 
centre. 

3.4  Variance  of  mean  slope 

The  derivation  of  the  variance  of  the  mean  water  height  given  in  section 
3.1  is  sufficiently  general  that  it  can  be  applied  with  only  minor  changes 
to  slopes.  By  analogy  with  equation  (16),  the  wave-number  spectrum  of  the 
mean  slope  of  the  illuminated  water  surface  is 


%Ck) 


2(k) 


2Ji  (kR) /kR  ^ 


(24) 


When  the  sea  is  fully  aroused,  the  wave-number  spectrum  is  given  by 
equation  (8).  A derivation,  similar  to  that  used  to  obtain  equation  (19), 
shows  that  the  variance  of  the  mean  surface  slope  relative  to  the  total 
slope  variance  is 


— — 

Sj^(0)/S(0)  = 

ln(ki  /ko) 

iko  ^ 

2 Ji  (kR)  /kR 

(25) 


The  dependence  of  this  ratio  on  the  dimensionless  parameters  koR  and  ki /ko 
is  shown  in  figure  4. 

It  can  be  seen  from  figure  4 that,  as  the  spot  shrinks  to  a point 
(koR^  0),  the  variance  of  the  mean  slope  of  the  area  Sj^(O)  approaches  the 

total  variance  of  the  water  surface  S(0) . This  asymptotic  behaviour  is 
necessary  for  wave  slope  variance,  as  it  is  for  wave  height  variance  and, 
consequently,  figure  4 resembles  figure  1 in  this  respect.  Furthermore, 
when  ko  R IT  and  the  diameter  of  the  spot  approaches  the  wavelength  of  the 
longest  significant  wave  present,  the  variance  of  the  mean  slope  falls 
rapidly  when  ki /.ko  = 1,  as  does  the  variance  of  mean  wave  height  in  figure  1. 

However,  the  parameter  ki /ko  is  much  more  significant  for  the  mean  slope 
variance  than  it  is  for  the  mean  height  variance.  In  other  words,  the 
mean  slope  variance  is  determined  by  the  whole  range  of  wave-numbers 
contributing  to  the  slope  spectrum  and  not  just  the  low  cut-off  wave-number 
ko . The  difference  is  due  to  the  functional  forms  of  the  wave  spectra  in 
the  equilibrium  range  given  in  equations  (7)  and  (8) . Whereas  the  height 
spectrum  is  clearly  dominated  by  small  wave-numbers  or  long  wavelengths, 
the  slope  spectrum  is  only  weakly  biased  in  that  way. 

In  any  application  requiring  the  detection  of  a laser  beam  reflected 
from  the  sea  surface,  the  strongly  directional  nature  of  the  reflection 
can  present  problems.  The  reflection  will  fill  a solid  angle  of  a size 
determined  by  the  roughness  of  the  water  surface.  The  central  direction 
of  this  solid  angle  is  determined  by  the  mean  slope  of  the  illuminated  area 
of  water.  If  both  the  incident  laser  beam  and  the  receiving  telescope  are 


ERL-0015-TR 


- 8 - 


vertical,  the  mean  slope  of  the  illuminated  spot  must  be  small  enough  for  the 
reflected  cone  to  include  the  vertical.  A small  mean  slope  variance  is 
therefore  necessary  but  it  is  not  sufficient  for  reliable  detection  of  the 
reflection  of  a laser  pulse  from  the  sea  surface. 

3. 5 Variance  of  residual  slope 

The  residual  slope  of  the  water  surface  at  the  point  t + within  the 
averaging  circle  is  defined  by 

- IrW  , C26) 

which  is  analogous  to  the  definition  of  residual  height  in  equation  (20). 

A similar  argument  shows  that  the  mean  autocorrelation  function  of 
residual  slope  is  given  by 

s'j^Cr)  = S(r)  - S^(r)  . (27) 


which  is  the  analogue  of  equation  (22).  The  ratio  of  the  residual  slope 
variance  to  the  total  slope  variance  ^C0)/S(0)  is  shown  in  figure  5 as  a 

function  of  the  dimensionless  parameters  koR  and  ki /ko . 

The  asymptotic  behaviour  of  all  the  curves  is  as  required,  namely  the 
residual  slope  variance  approaches  the  total  slope  variance  as  the  size  of 
the  illuminated  spot  increases.  The  range  of  wave-numbers,  over  which  the 
slope  spectrum  is  saturated  and  which  is  indicated  the  parameter  ki /ko , 
clearly  is  of  great  importance.  Consequently,  figure  5 is  not  sufficient 
to  deduce  the  reliability  of  detection  of  a laser  pulse  reflected  from  a 
water  surface.  This  question  is  examined  in  greater  detail  in  the  next 
section. 


4.  LASER  BEAM  REFLECTION  FROM  THE  SEA  SURFACE 

The  reflection  of  a laser  beam  from  the  sea  surface  is  influenced  by  the 
properties  of  the  illuminated  area  of  water.  These  properties  are  now  used  to 
derive  the  amplitude  and  variability  of  the  reflected  pulse,  as  well  as  its 
time  delay  and  time  spread. 

4.1  Average  amplitude  of  reflected  pulse 

The  reflection  of  a diverging  beam  from  an  airborne  laser  is  shown 
schematically  in  figure  6.  When  the  transmitted  radiant  flux  is  F^,  the 

irradiance  at  the  sea  surface  on  a plane  normal  to  the  incident  light  is 
approximately 

F cos^  (t> 

E(0)  = — , (28) 

nh"* 

where  h is  the  aircraft  height,  is  the  solid  angle  of  the  diverging  beam, 
and  0 is  the  nadir  angle.  It  is  assumed  that  the  irradiance  is  uniform 
within  the  beam  and  that  the  divergence  is  small. 

The  radiance  L(0),  reflected  from  the  sea  surface,  at  a zenith  angle  0, 
can  be  found  using  the  "scattering  cross  section"  derived  from  the  geometry 
of  reflection  from  a water  surface  by  Cox  and  Munk(ref.ll,  equation  (9)). 
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When  the  receiver  is  adjacent  to  the  transmitter  (as  shown  in  figure  6), 
their  expression  becomes 

= 0.25  p(o)  p(<^)  sec'^  (l>  , (29’ 

EC0) 

where  P(o)  is  the  reflection  coefficient  at  normal  incidence  and  p(0)  is 
the  probability  of  the  water  facet  having  a tilt  <l>  at  the  correct  azimith. 
The  probability  given  in  equation  (9)  reduces  to 

P(<^)  = (27ra^)"‘  exp(  -%  tan^0/a^)  , (30; 

when  the  rms  slopes  in  the  crosswind  directions  are  equal,  that  is 


a = a - a . 
c u 


This  simplification  introduces  only  small  errors  compared  with  the  effects 
being  described  in  the  present  theory. 

The  radiant  flux  dF^  entering  the  receiver  of  area  A from  the  water 

surface  element  within  the  viewing  solid  angle  dco  is  given  by 

dF^  = L(0)  (Ah"^  cos^  4>)  (h^  sec^  <l>  dw)  , (31' 

where  the  first  bracketed  term  is  the  solid  angle  subtended  by  the 
receiver  and  the  second  is  the  projected  area  of  the  sea  surface  element 
normal  to  the  viewing  direction  (see  figure  6) . Combining  equations 
(28)  to  (31)  yields  the  normalized  differential  reflectance 


O dFj  _ P(o)  A.  sec^^ 
F^  doj  “ 4 h^  2ffa^ 


tan^0 


which  must  be  integrated  to  obtain  the  total  received  flux. 

Equation  (32)  is  integrated  over  the  solid  angle  within  the  field  of 
view  of  the  receiver,  which  is  assumed  to  equal  the  solid  angle  of  the 
transmitted  beam  ^ . This  cone  has  an  angular  diameter  5 and  its  axis 
is  tilted  by  an  angle  0o  to  the  nadir,  as  illustrated  in  figures  6 and  7. 
The  sides  of  the  cone  nearest  to  and  furthest  from  the  nadir  are  inclined 
at  the  angles 

01  = I 5/2  - 00  I and  02  = 5/2  + 0o  . (32 

The  solid  angle  element  at  the  angle  0 within  the  field  of  view  is 


do}  = 2a(0)  0 d0  , 


where 


a(0)  = vr 


0<0i  , 


\ — — j ' 1 ^ 

arcos  (00^  + 0^ -5^ /4)/ (20o0)  0i<  0<02 
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For  tilt  angles  less  than  5/2  (figure  7(a)),  the  integration  is  between 
0 and  02 > while  for  larger  tilts  it  is  from  0i  to  02  (figure  7(b)).  The 
ratio  of  received  to  transmitted  flux  can  be  written  as 


p(o)  (A/4h^n) 

R(0o,S,a)  , 

(35) 

where 

/ 

■ sec^0  / 

tan^0'\ 

R(0o,5,0)  = j 

27ra^ 

2(j2  J 2a(0)0d0  , 

(36) 

is  the  relative  reflectance  of  the  sea  surface. 

Two  limiting  cases  of  the  relative  reflectance  can  be  evaluated  when 
the  beam  is  vertical  (0o  = 0,  a = tt,  tan  0 0)  . 

Firstly,  when  the  water  is  calm  (u-»-o) 

R(0,6,0)  = 1 . 

In  general,  therefore,  R(0o,5,a)  is  the  amplitude  of  the  reflection 
relative  to  that  from  a mirror.  Secondly,  when  the  water  is  rough 
(a»  6) 

R(0,5,a)  = S2 /27ra2  . (37) 

Equation  (35)  then  reduces  to 

F /F^  = p(0)  A / STTh^a^  , (38) 

r t 

which  is  a useful  result  for  practical  applications.  It  shows  the 
received  signal  to  be  independent  of  the  beam  divergence,  therefore,  the 
most  useful  quantity  to  evaluate  is  R(0o5,u)  / fl. 

Figure  8 shows  the  dependence  of  the  relative  signal  return 
R(0o,5,a)  / on  roll  angle  4>o  for  several  values  of  rms  slope  o and  beam 
divergence  6.  IVhen  the  water  is  rough  (0  — 0.2),  the  beam  divergence  has 
little  effect  and  the  retu:m  signal  is  relatively  insensitive  to  roll.  In 
calm  water,  however,  a small  roll  angle  can  produce  an  enormous  drop  in 
signal.  Under  these  conditions  the  beam  divergence  is  significant.  The 
approximate  value  of  the  relative  signal  return  determined  from  equation 
(36) , which  is  valid  for  0o  = 0 and  6 , is  shown  in  figure  8 as  a solid 
circle  on  the  left  axis.  The  agreement  between  these  points  and  the 
curves  for  6 = 10  mrad  provide  confirmation  of  the  accuracy  of  the 
numerical  integration  used  to  obtain  the  curves. 

The  range  of  seas  observed  by  Cox  and  Munk(ref.ll)  correspond  to  a 
variation  of  a from  0.09  to  0.28.  This  range  of  conditions  would  produce 
vertical  signal  returns  varying  from  2 to  20  sr~* , independent  of  the  beam 
divergence  below  100  mrad.  The  signal  at  15°  roll  angle  would  drop  by  a 
factor  between  about  2 and  100  over  this  range  of  sea  states. 

Figure  8 provides  only  the  average  amplitude  of  the  signal  return. 
Fluctuations  in  this  amplitude  due  to  the  passage  of  waves  are  equally 
important  and  are  examined  in  the  next  section. 
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4.2  Variability  in  amplitude  of  reflected  pulse 

In  the  previous  section,  the  calculation  of  the  received  pulse  amplitude 
assumed  the  average  slope  over  the  illuminated  area  to  be  zero.  Swell, 
however,  causes  fluctuations  in  the  average  slope  that  depend  on  the 
diameter  of  the  illuminated  spot  and  on  the  wave-number  spectrum  of  the 
sea  surface.  These  slope  fluctuations  induce  fluctuations  in  the 
amplitude  of  the  received  pulse,  because  of  changes  in  the  reflection 
geometry. 

When  the  transmitter  and  receiver  are  vertical  but  the  mean  slope  of  the 
illiMinated  surface  is  s,  the  reflection  geometry  is  equivalent  to  a 
horizontal  mean  sea  surface  illuminated  by  a beam  inclined  at  an  angle  0o , 
where  tan  0o  = s.  The  amplitude  of  a pulse  reflected  from  an  inclined 
surface  can  therefore  be  read  from  figure  8,  or  calculated  from  equation 
(36). 

When  a threshold  detector  is  used  to  detect  pulses  of  variable  amplitude, 
it  is  important  to  know  the  probability  that  the  pulse  amplitude  exceeds  the 
threshold.  It  can  be  seen  from  figure  8,  that  the  probability  of  the 
relative  signal  return  exceeding  R(^o>S><^)  / ^equals  the  probability  that 
s < tan  00 • Assuming  an  isotropic  gaussian  distribution,  the  probability 
of  a mean  slope  s at  any  azimuth  is 


pCs)  = s CT~^exp(-J2s"^  / a^)  , 


(39) 


where  a is  the  rms  mean  slope  (c.f.  equation  (30)  at  a particular  azimuth). 
Integration  and  inversion  yields 


00 


arctan 


a(ln  q”*)^ 


(40) 


where  q is  the  "failure"  probability,  that  the  mean  slope  exceeds  tan  0o 
and  the  pulse  amplitude  is  below  threshold. 

The  threshold  signal  retiorns  corresponding  to  different  failure 
probabilities  can  be  determined  from  equation  (36)  or  figure  8 with  the 
aid  of  equation  (40).  The  results  for  the  failure  probabilities  0.1%,  1%, 
and  10%  are  presented  in  figures  9,  10,  and  11  respectively.  It  can  be 
seen  that  the  threshold  must  be  set  consideraibly  lower  than  the  peak  signal 
return  if  reliable  detection  is  to  be  achieved  over  a wide  variety  of  sea 
conditions. 

The  threshold  signal  returns  in  the  absence  of  swell  and  the  absence  of 
aircraft  roll  are  given  by  the  curves  for  cr  = 0.  These  threshold  curves 
are  identical  for  all  failure  probabilities,  because  there  is  no  variation 
in  the  amplitude  of  the  reflected  pulse  when  the  mean  slope  over  the 
illxaninated  area  remains  zero.  Furthtsrmore , these  threshold  curves 

correspond  to  the  signal  returns  shown  in  f injure  8 when  0o  = 0,  that  is  in 
the  absence  of  aircraft  roll.  This  follows  from  the  fact  that  figure  8 
was  produced  with  the  assumption  of  a zero  mean  slope  (o  = 0) . 

Each  set  of  threshold  curves  for  a fixed  beam  divergence,  show  a fall  to 
lower  values  as  the  rms  mean  slope  increases.  This  is  because  the 
amplitude  of  the  reflected  pulse  in  the  vertical  direction  becomes  more 
variable  as  the  swell  increases  and  tilts  the  reflected  cone  further  away 
from  the  zenith.  The  reduction  in  the  threshold  is  greatest  when  the 
reflected  cone  has  a small  solid  angle,  that  is  when  the  rms  residual  slope 
within  the  illuminated  spot  is  small. 

Beam  divergence  is  significant  when  the  rms  residual  slope  is  small,  that 
is  when  the  water  surface  is  relatively  smooth  within  the  illuminated  spot. 
Under  these  conditions  and  in  the  absence  of  swell  (d  = 0) , a higher 
threshold  is  obtained  with  a small  beam  divergence.  This  is  because  smooth 
horizontal  water  reflects  a laser  pulse  into  a narrow  vertical  cone,  which 
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is  collected  efficiently  by  the  receiver  when  the  divergence  of  the 
transmitted  beam  (6)  is  smaller  than  the  divergence  of  the  reflected  cone 
(o).  When  8>a,  the  outer  parts  of  the  transmitted  beam  are  reflected 
away  from  the  receiver  and  are  thus  wasted.  A lower  threshold  is  then 
required  for  reliable  detection. 

However,  as  the  swell  grows  in  magnitude  (a  increases),  while  the  water 
surface  remains  smooth  (a  j>mall),  the  situation  is  reversed:  the  higher 
threshold  is  obtained  with  a large  beam  divergence.  The  collection 
efficiency  with  a large  beam  divergence  (5>  o)  is  still  small,  but  as  the 
directions  of  the  reflected  cones  from  small  patches  of  water  within  the 
illuminated  spot  are  wobbled  by  swell,  some  parts  of  the  beam  still  reach 
the  receiver.  With  a small  beam  divergence  the  reflected  cone  from  the 
small  illuminated  spot  has  a greater  probability  of  missing  the  receiver 
altogether  - hence  the  lower  threshold. 

The  nominated  failure  probability  significantly  affects  the  threshold 
only  when  a strong  swell  is  present  (o  large).  As  mentioned  above,  when 
0=0  the  threshold  curves  are  independent  of  the  failure  probability. 

As  an  example,  consider  the  case  when  both  the  rms  mean  slope  and  the  rms 
residual  slope  both  equal  0.1.  Figures  9,  10  and  11  reveal  that  the 

threshold  signal  returns  are  0.6,  1.7,  and  5 for  the  failure  probabilities 
of  0.1%,  1%  and  10%  respectively. 

These  curves  are  sufficient  to  predict  a satisfactory  threshold  signal 
return  of  the  spot  size  is  known  and  if  the  rms  mean  slope  over  that  spot 
and  the  rms  residual  slope  within  it  are  known.  Since  such  information 
is  not  always  available  additional  assumptions  must  be  made. 

When  the  sea  is  fully  aroused  and  the  slope  spectrtim  is  saturated  between 
the  lower  and  upper  cutoff  wave-numbers  ko  and  ki  , the  variability  of  the 
reflected  pulse  can  be  analysed  further.  The  relative  threshold  signal 
will  be  redefined  as 


Q(q,D)  = R(4>o  ,5  ,o)  /n  , 


(41) 


where  q is  the  failure  probability  and  D = 2R  is  the  diameter  of  the 
illuminated  spot.  This  definition  is  possible  because,  for  given  values 
of  ko  and  ki  , the  relative  threshold  signal  is  a function  of  q and  D only. 


as  is  now  shown.  The  solid  angle  of  the  transmitted  beam  is 

S2  = R 8V4  , (42) 

as  may  be  seen  from  figure  6.  The  angular  diameter  of  the  beam  is 

5 = D/h  sec  00  , (43) 

and  in  the  calculations  described  below  it  will  be  assumed  that  h = 500  m 
and  00  =0.  From  equation  (27)  it  follows  that 

a(D)2  = s'j^CO)  = B ln(ki/ko)  - a(D)^  , (44) 


where  the  total  slope  variance  is  obtained  by  integrating  equation  (8). 

The  value  B = 4.6  x 10“^  assigned  by  Phillips (ref .2)  and  Wu(ref.5)  is  used 
in  the  present  calculations.  Finally,  from  equation  (25), 
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5 (D)2  = Sj^CO) 


rk, 

— — 

J ko 

k"‘ 

2J,  (.kR)/kR 

-"dk  . 


(45) 


This  establishes  the  validity  of  the  functional  form  of  Q(q,D)  in 
equation  (41) . 

Figure  12  shows  the  threshold  relative  signal  return  Q(q,D)  as  a 
fimction  of  the  diameter  of  the  illuminated  spot  for  the  failure 
probabilities  q = 0.1%,  1%  and  10%.  Three  sets  of  curves  are  shown,  all 
for  slope  spectra  covering  one  decade  of  wave-numbers,  but  with  different 
lower  cutoff  wave-numbers:  ko  = 0.1,  1,  and  10  m~*  . These  curves  were 
computed  from  equations  (41)  and  (36)  expressed  in  terms  of  q,  D,  ko  and 
ki  with  the  aid  of  equations  (42)  to  (45).  However,  the  curves  can  be 
determined  graphically  from  figures  4 or  5 and  9,  10,  or  11.  From  given 
values  of  ko , ki  and  D can  be  calculated  the  dimensionless  parameters  koR 
and  ki /ko , which  yield  o and  a from  figure  4 or  5.  Then,  depending  on  the 
value  of  q,  figure  9,  10  or  11  will  yield  the  appropriate  threshold  signal. 

It  can  be  seen  from  figure  12  that,  when  the  diameter  of  the  spot  is 
large  enough,  all  the  curves  have  the  same  asymptote,  which  is  constant 
imtil  D exceeds  about  20  metres.  The  common  asymptote  is  due  to  the 
common  value  of  ki /ko  ( = 10),  and  hence  the  common  value  of  the  rms_total 

slope.  When  the  diameter  becomes  large  enough  the  rms  mean  slope  (cr) 

approaches  zero  and  the  threshold  signal  (see  figure  11)  for  a given  value 
of  o depends  on  the  beam  divergence  S only.  Figure  11  clearly  shows  that 
increasing  beam  divergence  reduces  the  threshold  return  when  a = 0 
(particularly  when  a is  small).  This  accounts  for  the  fall  in  the 

asymptotic  threshold  in  figure  12,  when  diameter  of  the  spot  exceeds  about 

20  m. 

Figure  12  shows  that  the  diameter  of  the  illuminated  spot  needed  to 
achieve  a given  threshold  return  depends  markedly  on  the  lower  cutoff  wave- 
niimber  ko  . More  specifically,  the  thresholds  begin  to  fall  in  every  case 
when  the  spot  diameter  becomes  less  than  the  longest  significant  wave 
present  (koD  < 2ir). 

The  acceptable  failure  probability  is  also  a significant  parameter. 

For  example,  when  the  slope  spectrum  is  saturated  over  the  range  1 to  10  m , 
a spot  diameter  of  1 m produces  a 10%  failure  level  with  a threshold  return 
of  10  sr”*  . In  order  to  improve  reliability  to  a failure  level  of  0.1%,  it 
is  necessary  either  to  reduce  the  threshold  by  an  order  of  magnitude  or  to 
increase  the  spot  diameter  by  a factor  of  about  2.5. 

The  range  of  wave-numbers  over  which  the  slope  spectrum  is  saturated  has 
a significant  influence  on  the  threshold  as  shown  in  figure  13.  All  the 
curves  correspond  to  a failure  probability  of  1%.  It  can  be  seen  that  the 
asymptotic  threshold  falls  as  the  saturation  range  (ki/ko)  becomes  larger. 

In  other  words,  as  the  sea  becomes  rougher  the  reflection  becomes  more 
diffuse  and  hence  less  intense  in  the  vertical  direction. 

The  fall  in  threshold  when  the  spot  diameter  becomes  less  than  the 
longest  significant  wave  present  (koD^  ZTT),  which  was  noted  in  figure  12,  is 
also  evident  in  figure  13.  However ,~the  fall  is  much  more  rapid  when  the 
slope  spectrum  is  saturated  over  only  a small  range.  When  the  saturated 
range  is  large,  due  to  the  presence  of  many  ripples,  the  threshold  is  not 
very  sensitive  to  the  size  of  the  illuminated  spot. 

The  above  results  are  presented  differently  in  figures  14  and  15.  By 
using  the  lower  and  upper  cutoff  wave-numbers,  k©  and  ki  , a.s  axes,  each 
point  in  the  plane  corresponds  to  a particular  sea  state.  The  small  lower 
cutoff  wave-numbers  to  the  left  of  the  graph  correspond  to  the  long  wave- 
length swell  of  the  open  ocean,  while  the  right  of  the  graph  represents  the 
small  waves  of  restricted  waters  such  as  rivers  and  bays.  On  the  vertical 
axis,  the  large  upper  cutoff  wave-numbers  towards  the  top  of  the  graph 
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represent  the  small  ripples  generated  in  a strong  wind,  while  the  lower 
part  of  the  graph  represents  smooth  waves  that  occur  in  calm  air.  A 
variety  of  sea  conditions  Cein  thearefore  be  represented  by  an  area  of  the 
graph. 

The  restrictions  on  the  illuminated  spot  diameter  D and  the  threshold 
signal  level  needed  to  achie;ve  failure  probabilities  less  than  1%  and  0.1% 
are  depicted  in  figures  14  £ind  15  respectively.  For  a given  threshold  and 
diameter,  the  failure  probability  is  below  the  required  figure  to  the  right 
of  the  corresponding  cuarve.  For  a threshold  of  10  sr“*  , reliable 
detection  is  achieved  only  vifhen  the  saturated  range  of  wave-numbers  ki  /k®  is 
less  than  30.  As  the  spot  diameter  is  increased,  reliable  detection  is 
achieved  in  the  presence  of  longer  swells.  The  effect  of  reducing  the 
threshold  (for  example  with  a more  sensitive  detector)  is  to  achieve 
reliable  detection  when  the  wave  spectriam  is  saturated  over  a greater  range 
of  wave-naambers . If  reliable  detection  is  required  over  a wide  range  of 

upper  cutoff  wave-numbers  and  with  lower  cutoff  wave -numbers  as  low  as 
0.1  m~^ , a diameter  between  1 and  10  m is  required  together  with  a threshold 
less  than  2 sr"‘  . 

A comparison  of  figures  14  and  15  shows  that  the  failure  probability  can 
be  reduced  from  1%  to  0.1%  vifith  an  increase  in  spot  diameter  or  a reduction 
in  the  threshold.  The  change  is  greatest  at  small  values  of  the  lower 
cutoff  wave-number,  that  is  in  the  presence  of  long  waves  when  the  reflected 
pulse  has  the  greatest  amplitude  fluctuations. 

4.3  Water  height  resolution 

The  time  spread  of  a laser  impulse  reflected  from  the  sea  surface  is 
determined  both  by  the  variation  of  water  height  within  the  illuminated  area 
and  by  the  variation  of  the  slant  height  travelled  by  the  beam.  This  is 
illustrated  in  figure  16.  The  accuracy  of  a water  height  measurement  is 

limited  by  the  sum  of  these  two  errors  and  hence  on  the  sea  state  and  the 
diameter  of  the  illuminated  spot. 

The  variance  of  the  residual  water  height,  which  can  be  derived  from 
equations  (19)  and  (22),  is  given  by 

J2B(ko  “"  - ki-^)  z'j^C0)/Z(0)  . (46) 

The  slant  height  error  of  a beam  at  a zenith  angle  <j>  from  an  aircraft 
at  height  h is 


h(sec  0 - 1)  . 


In  the  small  angle  approximation,  when  the  square  of  this  error  is 
averaged  over  an  area  of  diameter  D,  the  variance  of  slant  error 

D'’/192tf  , (47) 


is  obtained.  The  total  rms  height  error  is  therefore 


A h 


J2B(ko-'‘-  ki"?)  Z'j^C0)/Z(0) 


+ 


D'*/192h^ 


(48) 


The  upper  limit  of  this  error,  which  occurs  when  kj^ko,  is  shown  in 
figure  17  as  a function  of  Iko  and  D. 
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In  calculating  these  values,  it  is  assumed  that  the  aircraft  altitude  is 
500  m and  that  the  axis  of  the  transmitted  beam  is  vertical.  When  k©  is 
large  the  first  term  in  equation  (48)  is  negligible  and  the  height  error  is 
proportional  to  . In  this  limit,  only  small  waves  are  present  and  the 
slant  height  error  dominates.  When  k©  is  small,  the  first  term  in 
equation  (48)  dominates.  It  can  be  shown  that,  in  this  limit,  the  height 
error  is  proportional  to  D.  In  the  presence  of  strong  swell,  therefore, 
wave  height  errors  are  larger  than  slant  height  errors  except  for  very 
large  diameter  spots. 


5.  DISCUSSION  OF  APPLICATIONS 

An  airborne  system,  such  as  an  altimeter,  wave  height  profiler,  or  water 
depth  sounder,  that  relies  for  its  operation  on  the  detection  of  a laser  pulse 
reflected  from  the  sea  surface,  has  several  system  parameters  available  for 
optimization.  One  is  the  aircraft  operating  height,  but  this  is  often 
determined  by  other  considerations.  Another  is  the  peak  output  power  of  the 
laser  pulse,  but  this  may  be  limited  by  the  laser  design.  Yet  another  is 
vertical  stabilization  of  the  laser  beam  against  pitch  and  roll  variations  of 
the  aircraft.  Apart  from  the  latter,  the  only  system  parameter  available  for 
optimization  in  practice  may  be  the  diameter  of  the  illuminated  laser  spot  on 
the  sea  surface,  which  is  controlled  by  the  divergence  of  the  transmitted  laser 
beam. 

The  theory  of  laser  beam  reflection  developed  in  the  previous  sections  has 
been  concerned  with  the  influence  of  laser  spot  diameter  on  accuracy  and 
reliability  of  surface  detection.  Reliability  is  affected  by  fluctuations  in 
the  amplitude  of  the  received  pulse,  due  to  variations  in  the  mean  slope  and 
surface  roughness  of  the  illuminated  spot.  Accuracy  is  affected  by  pulse 
broadening  due  to  water  height  variations  within  the  illuminated  spot  and  by 
pulse  to  pulse  timing  variations  due  to  swell.  The  choice  of  a suitable  spot 
diameter  depends  on  both  the  application  and  the  range  of  sea  conditions 
expected. 

The  theory  described  in  this  report  has  been  developed  on  the  assumption  that 
the  wave  slope  spectrum  of  the  sea  is  saturated  between  two  cutoff  wave-numbers. 
Unfortuaately,  experimental  data  to  test  the  validity  of  this  assumption  is 
lacking.  Moreover,  the  lack  of  experimental  data  precludes  any  detailed 
predictions.  Cutoff  wave-numbers  can  be  deduced  from  slope  variance  data,  but 
vastly  different  results  are  obtained  depending  on  the  assumptions  employed,  as 
the  author  has  shown  elsewhere(ref .7) . Here,  therefore,  different  applications 
can  be  discussed  only  in  general  terms. 

A wave  height  profiler  must  illuminate  a spot  that  is  small  enough  to  profile 
the  waves  but  large  enough  for  reliable  detection  of  the  reflected  signal.  The 
optimum  diameter  will  depend  on  the  wave-number  of  the  longest  waves  present, 
the  height  resolution  required,  and  the  maximum  failure  probability  that  can  be 
tolerated.  Figures  14,  15  and  16  provide  the  information  necessary  to 
determine  the  optimum  diameter. 

An  altimeter,  required  to  measure  height  above  mean  sea  level,  must 
illimiinate  a spot  that  is  large  enough  for  reliable  detection  of  the  reflected 
signal,  large  enough  to  measure  mean  sea  level,  but  small  enough  to  contain 
slant  height  errors.  The  optimum  diameter  will  be  much  greater  than  that 
required  by  a wave  height  profiler.  Figures  1,  14,  15  and  16  provide  the 
necessary  information. 

A water  depth  sounder,  which  uses  mean  sesa  level  as  a reference  surface, 
requires  a spot  that  is  large  enough  for  reliable  detection  of  the  reflected 
signal  but  small  enough  for  adequate  accuracy  in  the  height  measurement.  The 
large  spot  suitable  for  an  altimeter  may  not  establish  the  mean  sea  level  with 
sufficient  precision  for  using  as  a reference  in  depth  sounding.  On  the  other 
hand,  the  small  spot  suitable  for  a high  resolution  wave  height  profiler  may 
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compromise  the  high  detection  reliability  needed  by  a water  depth  sounder. 
Consequently,  the  optimum  diameter  for  an  airborne  depth  sounder  may  lie  between 
the  other  two.  Again,  figures  14,  15  and  16  provide  the  relevant  information. 


6.  CONCLUSION 

The  reflection  of  a pulsed  laser  beam  from  a limited  area  at  the  sea  surface 
has  been  analysed  theoretically.  Temporal  and  amplitude  fluctuations  are  shown 
to  depend  on  the  wave-number  spectrum  of  the  sea  surface. 

Using  a two  parameter  characterization  of  the  wave-number  spectrum,  an 
expression  has  been  obtained  for  the  probability  of  the  reflected  signal  failing 
to  exceed  a given  threshold  in  the  detecting  system.  Results  calculated  from 
this  expression,  which  are  presented  graphically,  allow  the  optimum  diameter  of 
the  illuminated  spot  on  the  sea  surface  for  a given  application  to  be  determined. 

The  main  hindrance  to  immediate  application  of  these  results  is  the  lack  of 
experimental  measurements  of  the  two  parameters  on  which  the  theory  is  based. 
When  the  necessary  data  becomes  available,  the  results  can  be  used  in 
applications  of  lasers  for  wave  height  profiling,  altimetry,  and  water  depth 
sounding . 
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Figure  5.  Dependence  of  the  relative  variance  of  the 
residual  slope  on  the  dimensionless  radius 
of  the  spot 
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Figure  10.  Dependence  of  the  relative  threshold  signal 
return  on  the  rms  residual  slope  for  a 
failure  probability  of  1% 
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Figure  11.  Dependence  of  the  relative  threshold 

signal  return  on  the  rms  residual  slope 
for  a failure  probability  of  10% 
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Figure  12.  Dependence  of  the  relative  threshold  signal 
return  on  the  spot  diameter  for  three 
failure  probabilities 
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Figure  13.  Dependence  of  the  relative  threshold  signal  return 
on  the  spot  diameter  for  a failure  probability  of 
1%  and  a variety  of  sea  states 
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Figure  15.  Range  of  sea  states  for  which  given  spot 
diameters  and  relative  thresholds  achieve 
a failure  probability  below  0.1% 
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Figure  16.  Variation  o£  water  height  and  slant  height 
errors  over  illuminated  spot 
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